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I . Introduction 

The purpose of this study is to perform the following tasks: 

TASK A: Verify the basic mirror design parameters and performanc 

assessments of the NOAA-MSFC GOES x-ray telescope as 
provided to the principal investigator by MSFC technical 
personnel. 

TASK B: Perform computer ray trace analysis to establish spot 

size and point response function on-axis and off-axis 
at 5,10,16, and 20 arc-minutes by plotting the line 
spread function. 

TASK C: Determine resolution improvements and efficiency losses 

by shortening the hyperboloid length of the telescope 
to introduce intentional vignetting. 

TASK D: Generate a computer magnetic tape of the mirror radius 

as a function of position along the optical axis to be 
used in the mirror diamond turning activity. 

The general telescope system was assumed to be a paraboloid-hyperboloid in 
a Wolter Type 1 configuration. A brief discussion is given in section II 
of the equations which specify the telescope parameters and the resolution 
as a function of the collecting area using the empirical formulas of Van- 
Spreybroech and Chase. ^ Section III contains a discussion of the computer 
ray trace code used to determine the spot size and point response function 
for off-axis rays. In section IV, the measured resolution of the Goddard 
ATM x-ray telescope (S-056) is compared to the rms blur circle radius and 
the full width half maximum of the line spread function. An empirical 
scaling formula, Eq. 26, which transforms the rms blur circle radius into 
a more accurate measure of resolution, is introduced. The geometrical 
imaging properties of the proposed NOAA-MSFC GOES x-ray telescope are 
discussed in section V. The conclusions of this study are presented in 
section VI. Appendix A contains a discusuion of alternate mirror designs 
for the NOAA-MSFC GOES x-ray telescope. 
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1 1 . Defining and Analysis Equations of a WOI.TKR typo I X-ray Teles c ope 

In this section a brief discussion will he given of the equations 

that specify the mirror surface parameters of a WOLTFK 1 x-ray telescope 

and the resolution of the telescope as a function of the collecting area. 

2 

The notation is that of Mangus and Underwood. 

The defining equations of the paraboloid (l 1 ) and hyperboloid (II) 
surfaces are given by 


y P 2 - ■f(*x ? +'P) 

a 1 b 1 - 


(la) 


(lb) 


where the coordinate system is defined in Fig. 1. The constants a, b, 

c, and p are determined by requiring the paraboloid and hyperboloid to 

be confocal (see Fig. 1) and the glancing angle on both the P and 11 mirror 

to be © at the intersection point of the mirrors. Explicitly, 
m 

require the slope of P at intersection point (x v . ) to he 

p ,mm , • p ,nun 

tan 0 

m 


^ Y?> 'H "' 1 

i 






or 

Sp — \J ~f~(Xn Q 


( 2 ) 



If wo define the focal length, f, of llio telescope an the distance along 


tlio ray from the Intersection point of 1* and II mirrors to the focal 
point, then 

= fs,fir W = >W 

where it has been assumed that the slope of H at the intersection point 

is tan 3 Q and glancing angle of ray on H at the intersection point 
m 

is @ such that the ray will make angle 4 Q with the optical 
m w m 

axis. Evaluating x^ from (la) and using (3) gives 


Y — A *f CoS (3Sm) 

/x f? m i n 


Note further that if the ray reflected from P at intersection is extended 

to the origin, this virtual ray will make an angle 2 0 with both the 

f re*i m 

optical axis and the actual rav reflected 11 at intersection to the focal 

A 

point. Hence 


f 


ot' — 2.C 


where the last equality follows since OF is the distance between the two 
foci of the hyperboloid. One also knows from the general properties of 
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a hyperboloid 

2 L z - r^ 

C l + b ^ 

( 6 ) 

Hence, it remains to relate the constant a to f and & in order to define 

in 

completely the P and H mirror parameters. 

Since x,. ■ x , and y , ■ y - at the intersection point of 

H,max p,min p,min 'p,max v r 

the P and H mirrors, then Eq. lb at the intersection point can be rewritten as 


L^'fC 0 S ! 'Cae»,) - c] -f Sin (t0nj _ | 

o. 1 b‘ 

where Equ. 3 and 4 have been used. Equation 7 may be rewritten as 

•f Z [ *4 C4>sC^g»i)~ 1 l* 3 _ £ _ , 

T ^ ~ £1 - a* 

where Eqs 5-6 have been used. If the constant a is given by 

G. ~ ~ jj2, C^S (10m) " 1 'J 

then equation 7a is satisfied. The P and H mirror are given by 


(7) 


(7a) 


( 8 ) 


~f ~ /p,m,n ' h2n 

a - 

b l = c 1 -^ 

c - -F/a 
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and the focal length is defined by 


“fs**'* (tj&rn) “ Y p) nt/* “ /ri, mau) £ 


( 10 ) 


Referring th Fig. 1 it is clear that the axial length from the intersection 
of the P-H mirrors to the focal point is 




(ID 


where Eq. 4 has been used. 

The collecting area for axial rays incident upon the telescope is defined 
by 


A ““ ^ (YfrmcLx Y ^ 


(12) 


In the present application a specific value for y is determined by requiring 

p ,max 

the total length of the telescope, Jl to be given. That is. 


- 0 X 

^ T //x 
- (5K f)W< x x + 'f^J 


X 

y ?, **** 


(13) 


Also, the axial length of the P mirror is given by 




The length of the H mirror is determined by 


^ = X «, 




- X 


ht, m \r\ 


(14) 


( 15 ) 


ay 


where x^ is evaluated from the ray trace by evaluating where the r 

which is reflected from (x , x ) intersects the H mirror. The ray 

p,max p.raax 




•PWSj^ju iuw°'^. 
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trace code is discussed in section III. 

The empirical formulas of VanSpeybroech and Chase for the rms blur 
circle radius, (T^, are useful for evaluating alternative mirror designs. 
However, it is desirable to express <T Q explicitly in terms of the collecting 
area A. From Eq. 4 of Ref. 1, ^" D for a flat image surface is given by 


i "W’e A 

v d ~ 3,0 71 -tc»?e m Zp 


*2 


(16) 


where ^ is the off-axis field angle and 0*^ is measured in radians. The 
collecting area is approximately given by 


A * ** //>,„ 4 tun &n> 


(17a) 


or 


A - a V / S I n fafijia* 8" 


(17b) 


6 «. 


By comparing Eqs 16-17 , it is clear that^l6 i not an explicit function of A, 
since A depends on in Eq. 17b. In order to plot (7**^ versus A, it is necessary 

to obtain an expression for $ p in terms of A, and & n * From Eq. 17a 

we have 

0 - - — -r — - 

f My P ,m, n ta " e ' 


7 m 


(18) 


Eliminating JL p between Eq. 14 and 18 gives 


A 


yl r - f [<2 C0SteB m )-!j + 27) f fern Sm iH9~) 


( 19 ) 


which may be rearranged into a quadratic equation for the focal length f : 


-l( • 


f ~ ° (20 , 

which has the solutions 

[ - a _ ¥Ab**Vt*m)-i 1 ?' /x 
~ L T £71 "fan S ,rl J (2U 

a [2 6 »sm«m) ] 

The "+" sign in Eq. 21 corresponds to the physically valid solution. Therefore, 
Jtj is given by 

o - sir J./* 1 - flfjCatwW-Q V'*. 

T + 2 1 7 71iZ ” 9 m J ( 22 ) 

Equation 22 * represents the desired expression for £ expressed In terms 
of 9. j., A, and Usingj[ 22^0“ ^ may be evaluated from^l6 . Further 

discussion of Eqs. 16 and 22 will be given in Appendix A. 
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ill. Ray Tract? Co de* 

In this section, a brief discussion of the ray trace code used in this study 

is given. A detailed discussion of the ray trace equations as applied to a 

3 

Wolter Type 1 x-ray telescope has been given by Forman, et. al. The computer 
ray trace program can be broken down into the following parts: 


l.a. Define input constants : a,b,c,p,y , ,y _ ,.x jC„, 6. 

P jUlin p p »T1U!X II tQ 

b. Assign directions cosines to an Incident ray. It has been assumed 
that the incident ray is in the x-y plane and makes an angle d 
with respect to the x axis. The angle Q assumed values 0,1,2..., 
20 arc-mins. 

c. Set up a grid on the entrance pupil, which is an imaginary plane 
perpendicular to the optical axis located at x , such that 

p y IH3X 

each ray will pass through equal areas on the entrance pupil. For 
examt le, the polar coordinates R, Aon the entrance pupil are 
given by , _ , 

(p =■ O, *“/ 

where Ap> m 360/ (NPH0-1) . The radial coordinates are given by 
2 2 

where A r - (y - y , }/(NR0 - 1). When evaluating the 

p f IQdX pySQln 

rms blur circle radius from the ray trace data, it is sufficient 
to take NPHO ■ 140 and NRO ■ 30 for a total of 4,200 rays to be 
traced through the telescope at normal incidence. However, when 
the point spread function is being evaluated from the ray trace 
data, a total of 300,000 rays were traced through the x-ray tele- 
scope. As a result of the large number of rays traced for the 
point spread function, the angular variable P in the entrance 
pupil only assumed values from 0 to 180° with NPHO ■ 3,000 and 
NR) ■ 50. The symmetry of the telescope about the x-y plane is 
used to obtain the intercepts of the rays which would have passed 
through the entrance pupil for • I80°to 360°. 
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5 

f 


2. For each field Angle Q , separate rays are traced through each grid 

point (R, ) on the entrance pupil. The intercept (x^, y^) on the 

F mirror is evaluated. It is confirmed that the point (x ,y ) is on 

the physical P mirror surface of length X p* A vector fora of Snell's 

law is used to determine the direction cosines of the reflected ray 
from the P mirror. The intercept (x^y^) with the H mirror is 

computed and Snell's law for reflection is applied to evaluate the 
direction cosines of the ray reflected from the H mirror into the 
image space. The ray is accepted by the telescope if the point 
<XH,y,j) i* on the physical H mirror surface of length £ |j. The 

The ray intercept with the Gaussian ir.age plane is also computed. ^ 
The various summations required in the rms bins circle calculations 
are also evaluated. 

3. After completing the ray trace for all grid points at a given field 
angle, the rms blur circle radius is evaluated for a series of image 
planes by using the formulas of Ref. A. Al»o, the location of the 
optimum image surface, ie, loci of image points with minimum rms blur 
circle radius, is also computed. 

4. The point response function is evaluated by setting up a 40 x 40 
uniform grid on the imjge surface. For each field angle, the ray 
intercepts with the image plane are sorted into different image p.rid 
locations. The result is a distribution of rays over the image grid 
cells. The number of rays per image plane cell time the 

element of area A A (■ collecting area/total number of rays incident 
upon the telescope) is a measure of the point response function. 
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IV • Measured Resolution of the S-056 X-ray Telescope 

It has been observed by Werner’’ that the full width half maximum (FWHM) 

of the line spread function (LSF) is a more accurate measure of the observed 

resolution of an x-ray telescope than the rms blur circle radius. The 

experimental results which support this conclusion have been known for a 

number of years, but have not been generally recognized. It follows from 

3 

the point response function for S-056 x-ray telescope that the FWHM is 
approximately equal to 1 arc-sec at a field angle of 16 arc-min. However, 
from Fig. 2 the rms blur circle radius for S-056 as computed by the empirical 
formulas of VanSpeybroech and Chase, Eq. 16, is 8 arc-sec at a field angle 
of 16 arc-min. Whereas, the actual, measured resolution of the S-056 tele- 
scope for visible light is given by then's in Fig. 3. It follows by comparing 
Figs. 2 and 3 that the rms blur circle radius does not provide a reasonable 
measure of the resolution of an x-ray telescope. 

Since the rms blur circle radius is easier to evaluate than the point 
response function from ray trace data, it is desirable to establish a scaling 
technique between the rms blur circle radius and measured resolution. If one 
assumes the measured resolution, R, in arc-sec of an x-ray telescope is a 
linear function of the field angle , 0 , in arc-min which passes through the 
origin 

(23) 

then, fitting the data of Fig. 3 in a least squared sense yields the constant 

^ ■ 8.5375 x 10 arc-sec/arc-min. (23a) 

It is desirable to evaluate a scaling function a(Q) which will transform the 
rms blur circle radius into a meaningful measure of resolution. Assume 


S-056 X-PAY TELESCOPE 



FIELD A*1GLE ( arc-min ) 



FIFLD ANGLE ( arc-min. 

Fig. 3 
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R - >(9)^ 


(24) 


The two simplest forms for a( & ) are 

a(6) m <KQ (24a) 

a(O) - i/10 (24b) 

Evaluating o( by requiring c( & fit the (R/ 0*^) data in column 4 of Table 1 
in a least squared sense yields 

< K ■ 0.015635/ arc-rain (25a) 

Whereas, evaluating by requiring io fit the ((T D /R) data in column 5 of 
Table 1 in a least squared sense yields 

/ - 0.377207 / arc-min (25b) 

Table 2 displays *e<r D 

and * v ^ ) as a function of ^ . By comparing 
columns 2 and 3 of Table 2 with column 2 of Table 1, it is clear that column 3 
of Table 2 is a better fit to R. Hence, the empirical formula 

R = 0^/(0.377207 Q) (26) 

yields a better estimate of the measured resolution of the S-056 x-ray telescope 

than the rms blur circle radius data. From Fig. 3, it is apparent that the 

empirical formula, Eq. 26, offers a reasonable measure of the resolution of an 

x-ray telescope. We shall assume Eq. 26 is valid in general for x-ray telescopes 
will 

and apply Eq. 26 to the GOES x-ray telescope to establish the resolution from 

A 

the rms blur circle radius data. The point spread function for the GOES x-ray 
telescope is also evaluated in section V. 




TABLE 1 


1 

Q(<x* o*,,*) Q^fac-uc) 


.16875 


.33750 


.50625 


.6750 


.84375 


1.0125 


1.18125 


1.350 


1.51875 


1.68750 


.22884 


.67844 


1.34882 


2.23996 


3.35188 


4.68457 


6.23804 


8.01229 


10.00733 


12.22316 
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TABLE 2 

0 (OAC- mil) 


.00716 

.04243 

.12653 

.28017 

.52407 

.87892 

1.36544 

2.00435 

2.81642 


oO^p (gac - S 


3.82218 


30333 


.44965 

.59597 

.74229 

.88860 

1.03492 

1.18125 

1.32757 

1.47389 


1.62022 
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V . Geometrical Imagine Properties of the NOM-MSFC C O ES X-Ray Telesc ope 


Based on discussions with MSFC personnel, the mirror parameters were 
evaluated for a telescope defined by 

~f dy\. 

130 e^rc- mirt s. 

^ T =- Zoir\> — 76.2 cm. 

From Eq. 9, the P-H mirror parameters were found to be 


CL~ 

b ~ 

C - 


0 , ZKoV-VlS c '*- 
32. 83 l Z IH Cry >- 

2^^5 8707 cm. 

33 * 03 . Cr »• 


and from Eqs. 3-4,10-12 

✓ = 131.33 5 7 37 cm. 

J-”"\ q.tSlM'”- 

x l4a - 

* f ' = l°- 3 55 973 C".- 

yp. <*>»* ^ 

A = 25. 8179 66 c*,- 

aA = 8 . 6 >om- 3 *t 0~ 5 


Cm 


Yr^ 


(27) 


(28a) 


(24b) 


The length of P-H mirrors and,/ ~ are given by: 

F 


Xp - /#• 7/V0&3 cn >. 
/ = f. 30883 ¥C">- 

j = 65.285937 cm. 


(28c) 


The mirror parameters given by Eqs. 28 a »b,c have been evaluated independently 
by the principal investigator and are in agreement with those obtained by 
MSFC personnel. 
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A computer ray trace analysis has been performed on the telescope 
defined by Eqs. 28 to determine the rms blur circle radius as a function 
of the field angle for different image plane locations. Figure 4 presents 
the rms blur circle radius vs. the field angle for 7 image plane locations. 

The label " A z " in Fig. 4 specifies the displacement of the image plane from 
the focal point towards the telescope in cm. From the curves presented in Fig. 

4, it is clear that the depth of field is sharply defined and may present 
tight tolerances for positioning the image plane during assembly of the tele- 
scope. Figure 5 displays the resolution as predicted by the empirical formula, 

Eq. 26, for the GOES x-ray telescope. Based on these results, one may expect 
an actual resolution for the GOES x-ray telescope to be less than 2.5 arc-sec 
at a field angle of 20 arc-min. In order to confirm these results, the 
meridional (x-y plane) line spread function (LSF) and the sagittal 

LSF have been evaluated from the point spread function as computed from the 

arc-min 

ray trace data for field angles Q ■ 5,10,16,20^. The meridional LSF's are 

arc-min 

given in Figs. 6a,b,c,d for £) ■ 5, 10, 16, 20^, respectively. The FWHM's of the 
LSF's given in Figs. 6a-d are displayed in Fig. 5 as the symbols. A 

straight line which approximates the FWHM of the meridional LSF has also been 
given. (The plate factor used to express the FWHM in angular units is 
(180*360/^" f) ■ 3,123.314 arc-sec/radian where f ■ 66.04 cm for GOES.) It 
may be concluded from Fig. 5 that the empirical formula, Eq. 26, and the FWHM of the 
meridional LSF agree reasonably well at large field angles and differ by a 
factor of 2 for small field angles. Whereas, the rms blur circle radius and 
the FWHM of the LSF differ by almost an order of magnitude. 

The sagittal LSF's for a field angle of 5,10,16,20 arc-min are given in 
Figs. 7a,b,c,d. The FWHM's for Figs. 7a-d are plotted in Fig. 5 as the "+" 
symbols. From these results. Figures 5-7, it follows that the resolution of 
the GOES x-ray telescope should be less than 3 arc-' ee at a field angle of .• 
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20 arc-min and less than 1 arc-sec for a field angle of 5 arc-min. It is 
recognised that the accuracy with which the empirical formula, Eq. 26, predicts 
the resolution of an x-ray telescope can he improved by correlating the observed 
resolution with the rms blur circle radius data for additional x-ray telescopes* 
This has not been done due to time constraints under the present study. However 
it should be noted from Fig. 5 that the empirical formula * Eq. 26* agrees 
reasonably well with the FWHM results. 

The discrepancy between the resolution predicted by the rms blur circle 
radius and the FWHM from the LSF may be understood by examining the point spread 
function data in more detail. It has been found that for a field angle of 16 
arc-min only 30% of the incident rays intercepts the image plane within a 
region centered on the central peak in the point spread function of a diameter 
equal to the FWHM. The majority of the rays are incident upon the image plane 
at considerably larger distances from the central peak in LSF. Thus* the rms 
blur circle radius is larger than the FWHM. However, the high energy density 
located within the central peak of the LSF dominates the spread out backgound 
illumination at the detector in producing an image. Thus, the FWHM of the LSF 
is a more accurate measure of resolution of an x-ray telescope than the rms 
blur circle radius. 

2 

It has been observed by Mangus and Underwood that the off axis resolution 
as measured by the rms blur circle radius of an x-ray telescope may be improved 
by shortening the length of the H mirror. Figures 8-13 present the rms blur 
circle radius vs. the field angle for / H - 9.0,8.75, 8.5, 8.25, 8.0, 7.5, 
respectively. By comparing Figs. 2 and 13, it is seen that the resolution 
may be improved by approximately 10% by shortening the length of the H mirror 
to 7.5 cm. and defocusing the image plane by +0.2 mm. However, the collecting 
area of the telescope is reduced. From Fig. 14, it is seen that when j£ ^ 


is 
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Fig. 10 
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r educed to 7.5 cm the collecting area is reduced by approximately 25% to 

2 

19 cm . Since the reaolution of the proposed NOAA-MSFC GOES x-ray tele- 
scope defined by Eqs. 28 is less than 3 arc-sec at a field angle of 20 
arc-min, it should not be necessary to introduce vignetting of the H mirror* 
to defocuse the image plane* or to use a curve image surface in order to 
meet the performance objectives of the GOES x-ray telescope. 

In order to complete the discussion of the geometrical imaging properties 
of the GOES x-ray telescope with JL 9.3089 cm, the rms blur circle radius 
as a function of the field angle on a flat and curved image surface is 
given in Fig. 15a. Figure 15b presents the rms blur circle rauius as a function 
of A z of the optimum image surface for the f.eld angles 4*8*12*16*20 arc- 
min. Figure 15c displays A s of the optimum image surface vs the field angle. 
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VI . Conclusion 


It has been shown (Fig. 5) that a resolution better than 3 arc-sec over 
a + 20 arc-min field of view should be obtained for the NOM-MSFC GOES x-ray 
telescope defined by Eqs. 27-28. It has also been shown (Fig. 13) that a 
10Z improvement in the rms blur circle radius results when is shorten 
to 7.5 cm and when the image plane is defocused by 0.2mm towards the tele- 
scope. However, the collecting area of the telescope is reduced by 25Z 
when is shorten to 7.5 cm (Fig. 14). The effect of the intentional 
vignetting of the H mirror on the FWHM of the LSF has not been explored. 

The shape of the optimum curved image surface based on the rms blur circle 
analysis has been evaluated (Fig. 15c). 

TASK A-C have been fully addressed and the results are presented in 
sections II-V of this report. It is the conclusion of this study that 
the proposed NOAA-MSFC GOES x-ray telescope defined by Eqs. 27-28 should 
meet the performance objectives for the resolution and collecting area 
which have been setforth by MSFC technical personnel for this project. 

In order to complete TASK D, it is necessary for MSFC technical personnel 
to provide the format information required to generate the computer magnetic 
tape of the mirror radii as a function of the position along the optical 
axis for use in the diamond turning activity. Since the format information 
for diamond turning activities is not currently available from MSFC, the 
principal investigator proposes to complete TASK D after the submission of 
the final report, as soon as the format information is delivered to the 


principal investigator. 
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APPENDIX A: Alternate Mirror Designs for the GOES X-ray Telescope 

An analysis of the empirical formulas of VanSpeybroeck and Chase* for the 

resolution of a telescope on a flat image plane as a function of the collecting 

area is given in this section. Specifically, the rms blur circle radius, CT^, 

(Eq. 16 where Eq. 22 is used for JL ) has been plotted as a function of the 

collecting area for maximum glancing angles , O = 110,120,130,140,150 arc-min. 

m 

Figures A-la,b,c,d are plots of & vs the collecting area for field angles 

5,10,16,20 arc-min. In Fig. A-lc for & « 16 arc-min, it is seen that a rms 

blur circle radius of 20 arc-sec may be achieved at Q - 130 arc-min and 

m 

2 

a collecting area of 19 cm , which is consistant with the results of the 

extensive ray trace described in section V. Minor improvements may be expected 

over those contained in Figs. A-l when defocusing and vignetting of the H 

mirror are performed. For convenience, the focal length f from Eq. 21 as a 

function of A, 0 , and is presented in Fig. A-2. In order to translate 
m T 

specific values of A and 0 from Figs. A-la-d into corresponding focal length, 

m 

one uses Fig. A-2. 
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